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Abstract 

- - - Track detectors in high energy physics experiments require an accurate de- 

^ . termination of a large number of aUgnment parameters. A general method 

O ' has been developed, which allows the determination of up to several thousand 

(N . alignment parameters in a simultaneous linear least squares fit of an arbitrary 

bp. number of tracks. The sensitivity of the method is demonstrated in an example 

of the simultaneous alignment of a 56-plane drift chamber and a 2-plane sili- 
con tracker. About 1400 alignment parameters are determined in a fit of about 
fifty thousand tracks. 
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1. ALIGNMENT PROBLEMS OF TRACK DETECTORS 

Alignment problems for large detectors, especially track detectors, in particle physics often require the 

determination of a large number of alignment parameters, typically of the order of 100 to 1000. Align- 

5o \ ment parameters for example define the accurate space coordinates and orientation of detector com- 

O ■ ponents. In addition drift chambers require the determination of the drift velocity Udrift and of local 

^ , variations of the drift velocity, and of the Lorentz angle. Figure [T] shows an example of residuals in a 

' silicon vertex track detector before and after an alignment. 

X ' 

D \ In the alignment usually special alignment measurements are combined with data of particle reac- 

Q^' tion, typically tracks from physics interactions and from cosmics. In this paper the alignment parameters 

O . are called global parameters; they contribute to all the data. Usually these parameters are corrections to 

default values and are thus relatively small, and the value zero is the initial approximation. Parameter 
sets valid for a single track are called local parameters, examples are track slopes and curvatures. 

^ ' An approximate alignment method is to perform least squares fits on the data ignoring initially 

■ the alignment parameters. The deviations between the fitted and measured data, the residuals, are then 

used to determine the alignment parameters afterwards. Single fits depend only on the small number of 
local parameters like slope or curvature of the track, and are easy to perform. This approximate method 
however is not a statistically correct method, because the (local) fits depend on a wrong model (namely: 
ignoring the global parameters), and the result are biased. In addition the results may be unstable due 
to detector inefficiencies. The adjustment of parameters based on the observed (biased) residuals will 
then result in biased alignment parameters. If these alignment parameters are applied as corrections in 
repeated fits, the remaining residuals will be reduced, as desired. However, the fitted parameters will still 
be biased, even more, if this approximate procedure is applied iteratively. The simple residual method 
is also unable to determine correlated parameters, for example a geometrical shift of a detector position 
and a shift in drift velocity in a drift chamber. 

A better method is an overall fit, with all the (global) alignment parameters and local parameters, 
perhaps from thousands or millions of events, determined simultaneously. The most important general 
method for the determination of parameters from measured data is the linear least squares model. It is 
usually stable, fast and accurate, requires no initial values of parameters and is able to take into account 
all the correlations between different parameters, even if there are many parameters. There is of course 
a practical limit in the number of parameters in simultaneous fits, because of space and time limitations. 
In an alignment fit the size of the vectors and matrices in the least squares solution will be large, of the 
order of 10^ or 10^. With standard methods such a solution is impossible. 
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Fig. 1: Residuals between tracks and hits of a silicon tracker are shown as a function of the 
azimuthal angle (vertical axis, in degrees) before and after an alignment. The residual on the 
horizontal axis is given in units of 10 nm. 



In this paper a special method of solution for this kind of problems is derived. Due to the special 
structure of the matrices with one set of global parameters and many sets of local parameters the problem 
can in fact be reduced to a solvable size, without making any approximations. The rather general method 
is realized in a program called Millepede [Q]. 

2. STANDARD LEAST SQUARES 

Linear models. In a linear model the measured quantity has an expectation value which is a linear 
combination of the parameters a with fixed factors (e.g. from geometry) which are combined to a vectore 
dk- The difference between the measured value and the expectation value is the residuum r^: 

Zk = -dk + rk ■ (1) 

The measured data are assumed to be independent; the result of a measurement does not depend on 
any other measurement. Thus the covariance matrix of the measured data is a diagonal matrix, which 
simplifies the treatment. The accuracy of the measurement Zk is given by the variance a'j. or the standard 
deviation fx^. It is assumed, that the precision is at least approximately known. 

The solution of the method of least squares is determined by the minimum of the weighted sum of 
the square of the residuals 



(2) 



where the sum (index k) is over all measurements. This sum is a quadratic function of the parameters a, 
and the solution is obtained from a set of Unear equations (so-called normal equations of least squares), 
which can be written in matrix form 

Ca = h with C = uj^ (d^ • d^) and b = ''^^Wf^Zkdi^ (3) 

k k 
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with a symmetric n-hy-n matrix C and a right-hand-side n-vector 6; the elements of C and b are sums 
with independent contributions from each measurement (because the measured data are independent). 
There is no Umitation on the number of measured data; it is not necessary to keep all input data in 
memory, because the complete information from a single measurement is contained in the contributions 
to C and b, which have a fixed size. The solution vector is a = C^b with the inverse of the matrix C, 
which, by error propagation, is equal to the covariance matrix of the parameters a. 

What are the limits for the number n of parameters, if the matrix equation is solved by standard 
methods? The symmetry property of the matrix C may be used to store only the (n^ + n)/2 elements on 
and above the diagonal, and also to reduce the computer time to almost half the time for a non-symmetric 
matrix. Matrix inversion is an explicit process and requires a CPU time proportional to for a n-by-n 
matrix. For a 1000-by-lOOO symmetric matrix in double precision (about 4 Mbyte) the CPU time is a few 
minutes on a standard PC. The time (oc n^) and space (oc n^) requirements imply at present a practical 
limit of the size of the matrix (and number of parameters) between n = I 000 and n = 10 000. 

The accuracy of the computation of the inverse matrix is a critical point, and often computation 
in double precision is required, especially for large matrices. The accuracy is generally low for a large 
condition number of the matrix (the condition number is defined as the ratio of the largest to the smallest 
eigenvalue of the matrix). From the interpretation of the inverse matrix as covariance matrix it is clear, 
that the condition number will be large if the correlation between parameters is high. The matrix is 
singular (one or more eigenvalues are zero) in case of a complete correlation between parameters. In 
this case and also in cases of missing data (e.g. dead channels in an alignment fit) a standard matrix 
inversion program will fail. A singular value decomposition (SVD) would allow to recognize the cases 
of zero or small eigenvalues, which would destroy the full solution. The standard method in case of zero 
or very small eigenvalues is to set their inverse to zero; this means essentially to remove parameters or 
linear combinations of parameters from the solution (i.e. the resulting parameter values are zero); if the 
parameters are in fact corrections this is clearly acceptable. In the program Millepede iQ] a similar, but 
simpler and faster strategy is used to avoid the effects of a bad condition number. This strategy also 
works in cases, where due to missing input data|] one or more parameters are undefined; the parameters 
will become zero in this case. The strategy allows in addition to remove certain parameters from the least 
squares fit without the necessity to reorganize the numbering. 

Solution by partitioning. The special structure of the matrix C to be inverted may allow to move the 
limit of solvable problems further (see also [p). Below the matrix C is partitioned into submatrices, and 
the vectors a and b are partitioned into two subvectors; then the matrix equation becomes 
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where the submatrix Cu is a p-hy-p square matrix and the submatrix C22 is a q-hy-q square matrix, 
with p + q = n. If the sub-vector ai would not exist, the solution for the sub-vector a2 would be defined 
by the matrix equation 

C22 0-2 — b2 with solution 02 = C22 ^2 , (5) 

where the star indicates the special character of this solution. Only the q-hy-q sub-matrix C22 has to be 
inverted. The solution 03 of course differs from the final solution 02 and is called the local solution. 

Now, having inverted the sub-matrix C22, the submatrix of the complete inverse matrix C corre- 
ponding to the upper left part Cu, called B, is obtained by the formula 

B = [Cu — C12C22C12) , (6) 

'Large detectors may have dead channels, which means missing data. 
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which requires in addition to some matrix products the calculation of the inverse of a p-hy-p sub-matrix. 
With this matrix B the solution of the whole matrix equation can written in the form 
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which can be used to obtain the complete result (ai and a2). Note that to solve the matrix equation two 
smaller matrices, a symmetric p-hy-p and a symmetric q-hy-q matrix, have to be inverted. The sub-vector 
ai is obtained by the product 



ai = Bbi 



BC12C22 ^2 



B{hi-Ci2a*2) 



(8) 



which is simplified by using the special local solution 03 = C22 ^2- If the interest is the determination of 
this sub-vector ai only, while the sub-vector 02 is not needed, then only the equation (^) has to be solved. 
Some computer time can be saved by this method, especially if the matrix C22 is easily calculated or 
already known before. In special cases the method can be applied repeatedly, and therefore problems 
with a huge number of parameters become solvable. The method is not an approximation, but is exact 
and it takes into account all the correlations. This type of application is discussed in more detail in the 
next chapter. 

It seems that this possibility of removing unnecessary parameters, but still getting the correct 
solution for the remaining parameters in least squares problems was already known and has been used in 
the nineteenth century; one example is 

Constraints. If there are explicit relations between the parameters of a least squares fit, then they should 
be explicitly taken into account in the solution. The case of a single linear constraint between parameters 
of the form • a = /q is assumed here. The Lagrange multiplier method is the standard method to 
include constraints in a least squares solution; one additional vector A is introduced, with one element 
per constraint. The equation (|3|) is modified by the additional parameters to the form 
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where the matrix on the left hand side is still symmetric. Linear least squares problems with linear con- 
straints can be solved direcly, without iterations and without the need for initial values of the parameters. 



3. LEAST SQUARES WITH GLOBAL AND LOCAL PARAMETERS 

Local parameters. Local parameters are denoted by a^ , j = 1, 2, ... i^. A set of measurements z is 
considered. A single value z from the set is ideally given by a linear expression 2; = ai • (5i + a2 ■ ^2 + 
. . .ay • by, depending on the parameters and on (known) constant factors bj. More complete, the 
dependence can be written as 

V 

Zk = ai- bik + a2 ■ b2k + ■ ■ -a^ ■ buk = aj ■ bjk , (10) 

i=i 

with an additional index k indicating the k-th measured value. 

As an example the measurement of a straight particle track is considered, where each is a 
coordinate measurement at plane k. The local parameters are intercept (ai) and slope (02) of the particle 
track. If the distance of the plane is denoted by Sk, the linear track model can be written as 

Zk = ai ■ 1 + a2 ■ Sk . 
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From a least squares fit of the model to the set of measured z-values the local parameters intercept 
and slope can be determined. The formalism of least squares allows to take into account the different 
precision of the different measurements, by assigning a weight to each measurement ofwk = l/o"^. The 
least squares solution, minimizing the weighted sum of squared residuals, is defined by the so-called 
normal equations, which can be written in terms of matrices 

Ta = (3, solved by ct = . (11) 

The elements of the matrix T and the vector (3 are formed by sums 

Tij = '^Wk- 6ik5jk Pj = '^'Wk- ZkSjk ■ (12) 

k k 

These formulas are valid in the case of uncorrelated z-measurements (uncorrelated measurements are 
assumed thoughout this paper). The inverse matrix T^^ is the covariance matrix of the parameter vector 
ex. 

Global parameters. Now global parameters are considered, which contribute to all the measurements. 
They are assumed to be alignment parameters here, although the formalism is general and can also be 
applied to other problems. Aligment parameters are defined to represent corrections to ideal (design) 
values; this parametrization adds further terms to the equations (p^ 

n V 

z = ai ■ di + a2 ■ d2 + . ■ - an ■ dn + ai ■ 6i + 02 ■ 6 2 + ■ ■ ■ a„ ■ 6^ = aj ■ dj + ^ aj ■ 6j . (13) 

^ V ' ^ V ^ ._ ._ 

global parameters local parameters *~ J—'- 

Usually only few terms containing global parameters are nonzero for a single z measurement. In the 



following it is assumed that, as written in equation (13), the dependence on the global parameters is lin 



ear. A non-linear relationship for the global parameters would require iterations, starting with reasonable 
initial values for the parameters, and assuming a linearized expression in each iteration. 

The simultaneous fit of global and local parameters. In the following it is assumed that there is a set 
of partial measurements. Each partial measurement, with index i, depends on local parameters cti, 
and all of them depend on the global parameters. In a simultaneous fit of all global parameters plus local 
parameters from N subsets of the data there are in total {n + N -u) parameters, and the standard solution 
requires the solution of (n + • z^) equations with a computation proportional to (n + • z^)^. Below it 
is shown, that the problem can be reduced to n equations with a computation proportional n^. 

Generalizing the formalism of equations ( |To| ) to the complete measurement, starting from equa- 



tions of the type of (|13|), one obtains a system of normal equations with large dimensions, as is shown 
in equation (|l^). The matrix on the left side of equation ( |T4l ) has, from each partial measurement, three 
types of contributions. The first part is a contribution of a symmetric matrix Cj, of dimension n (number 
of global parameters). All the matrices Ci are added up in the upper left comer of the big matrix of the 



normal equations. The second contribution is the symmetric matrix Fj (compare equation (12)), which 
gives a contribution to the diagonal of the big matrix and is depending only on the i-th partial measure- 
ment. The third contribution is a rectangular matrix Gi, with a row number of n (global) and a column 
number of v (local). 

There are two contributions to the vector of the normal equations, bj for the global and /3j for the 
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local parameters. The complete matrix equation is given by 
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In this matrix equation the matrices Cj, F j, Gi and the vectors hi and /3j contain contributions from the i- 
th partial measurement. Ignoring the global parameters one could solve the normal equations ViCx* = (3^ 
for each partial measurement separately by 



(15) 



The complete system of normal equations has a special structure, with many vanishing sub-matrices. 
The only connection between the local parameters of different partial measurements is given by the 
sub-matrices Gi und Cj, 

The aim of the fit is solely to determine the global parameters; final best parameters of the local 



parameters are not needed. The matrix of equation (14) is written in a partitioned form. The general 
solution can also be written in partitioned form. Many of the sub-matrices of the huge matrix in equation 
([l4|) are zero and this has the effect, that the formulas for the sub-matrices of the inverse matrix are very 
simple. By this procedure the normal equations 



(16) 

are obtained, which only contain the global parameters, with a modified matrix C' and a vector b', 

C' = Y,C^-Y, G.TT^Gj 5' = ^ - ^ . (17) 



This set of normal equations contains explicitly only the global parameters; implicitly it contains, through 
the correction matrices, the complete information from the local parameters, influencing the fit of the 



global parameters. The parentheses in equation ( |17D represents the solution for the local parameters, ig- 
noring the global parameters. The solution a = C"^ h' represents the solution vector a with covariance 
matrix C'^. The solution is direct, no iterations are required. Iterations may be necessary for other 
reasons, namely 

• the equations depend non-linearly on the global parameters; the equations have to be linearized; 

• the data contain outliers which have to be removed in a sequence of cuts, becoming narrower 
during the iteration; 

• the accuracy of the data is not known before, and has to be determined from the data (after the 
alignment). 
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Application. A very important point in problems discussed above is the definition of the set of global 
parameters. It is essential to define them in a way which is sensitive to the data used to determine them. 
Computing the inverse of a huge (symmetric) matrix can be a delicate numerical problem. The matrix 
will be singular or close to a singular matrix, if one or more of the global parameters are not really 
determined (and constrained) by the data, and this could destroy the whole determination. This will 
happen in case of a strong correlation between different global parameters. 

It is wise to start with a few global parameters and to observe the correlations between them, 
before adding more global parameters. For example if too many global parameters are used in a detector 
alignment, the whole detector may shift and rotate freely in the fit. In those cases one could either fix 
certain global parameters or one has to introduce constraints (e.g. average displacement and rotation 
equal to zero). 

Experience with the method shows that it is essential to make a simultaneous alignment of all track 
detector components, which are also used for the measurement of tracks, while independent internal 
alignment of the single track detector components may not lead to a good overall result. It is also 
important to use, if possible, different data sets simultaneously. For example in a track detector alignment 
high-momentum cosmic muons may be useful, which traverse the whole detector and thus relate (and 
constrain) different parts of the detector. But this should be done simultaneously with tracks from physics 
events, because these events finally require a good alignment. 

The method described here has been realized in the program Millepede [^, written in Fortran??. 
It provides a set of subroutines for the mathematical methods and allows to adapt the method to a par- 
ticular problem. The program allows to adjust the dimension of vectors and matrices via PARAMETER- 
statements and to introduce linear constraints. Throughout the program use is made of the sparse charac- 
ter of the vectors and matrices in order to reduce the execution time. The program includes a method for 
outlier rejection, which in practice may be essential in alignment problems especially if some detector 
components have a bad initial alignment. This feature requires to iterate, realized with an intermediate 
data file, which is written during the first iteration and read during the other iterations. 

Sometimes the model underlying the alignment is nonlinear and also constraints may be nonlinear. 
The standard method to treat these problems is linearization: the nonlinear equation is replaced by a linear 
equation for the correction of a parameter (Taylor expansion); this requires a good approximate value of 
the parameter. 

4. ALIGNMENT OF THE CENTRAL TRACK DETECTORS IN THE HI EXPERIMENT 

The example described is the alignment in the rc^-plane perpendicular to the beam line of a 56-plane 
drift chamber and a 2-plane silicon vertex detector in the HI detector at HERA iQ]. Both are cylindrical 
detectors; the drift chamber has a length (z) of about 2 m, and extending from 20.3 cm to 84.4 cm in 
radius r. The silicon vertex detector [^] has two planes around the beampipe made of 0.15 mm Alu- 
minium and 0.9 mm carbon fiber. The drift chamber and the silicon vertex detector have a rip resolution 
of about 150 /im and 15 /xm, respectively. These detectors are interspersed by additional chambers for 
the z-measurement. 

The central jet chambers CJCl and CJC2 ([Q]) with an active length of 2200 mm and an outer 
radius of 844 mm have in total 2640 anode sense wires parallel to the beam with two adjacent cathode 
planes (wires) shaping the drift field. A jet chamber cell extends azimuthally from the sense wire plane to 
both adjacent cathode wire planes, and radially extends over the full radial span of CJCl and CJC2 with 
no further subdivision. The jet cells are tilted by about 30° such that in the presence of the magnetic field 
(1.1 Tesla) the ionization electrons drift approximately perpendicular to high momentum vertex tracks. 

The silicon vertex detector CST consists of two cylindrical layers of double sided, double metal 
silicon sensors read out by a custom designed analog pipeline chip. The two layers of the CST are formed 
from 12 and 20 faces at radii of 5.?5 cm and 9.?5 cm, respectively. One face or 'ladder' consists of six 
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row. 


number 


parameter 


a 


unit 


1 


2 


/S.X 


1 


fini 


2 


2 


Ax/Azr 


2 




3 


2 


Ay 


1 


fim 


4 


2 


Ay/Azr 


2 


/im 


5 


2 


Atp 


10 


fimd 


6 


2 


Atf/AZr 


10 


/irad 


7 


2 


AoLor 


100 


/irad 


8 


2 


AVtinii— /VdTift 


10"^ 




9 


2 


Avdrift+ / Vdrift 


IQ-^ 




10 


2 


ATb X Wdrift 


< 1 


fim 


11 


2 


wire staggering in wire plane 


few 


fim 


12 


2 


wire staggering perp wire plane 


few 


/im 


13 


2 


sagging in wire plane 


few 


fim 


14 


2 


sagging perp. wire plane 


few 




15 


180 


Andrift / f drift per cell half 


few 


10-* 


16 


112 


Audrift/fdrift per layer half 


few 


10-* 


17 


330 


ATo X Udrift per group 


10 


fim 


18 


56 


wire position in driftdir. per layer 


10 


fim 


19 


56 


ATo X Wdrift per layer 


10 


fim 


20 


56 


Axw per layer 


few 10 


fim 


21 


112 


Andrift/^^drift for Ie/50 UlA 


few 10"* 




22 


90 


Aiidrift/i'dnft per layer 


few lO"-* 




23 


90 


Ayw per layer 


few 10 


fim 


24 


90 


(Ayw)^ per layer^ 


few 10 


fim 


25 


64 


A in ladder 


few 


fim 


26 


64 


A perp. ladder 


few 


fim 


27 


64 


rel. A in ladder (zr) 


few 


fim 


28 


64 


rel. A perp. ladder (zr) 


10 


fim 


29 


64 


rel. A perp. ladder (ip) 


few 


fim 



Table 1: Alignment parameters determined in 
the Millepede fit. The Millepede accuracy is 
given by one standard deviation a. The differ- 
ent parameters are discussed in the text. 



silicon sensors and aluminium nitride hybrids at each end. A double layer of carbon fiber strips with a 
total thickness of 700 /lim and a height of 4.4 mm is glued to the edges. The position of the ladders in a 
layer are shifted tangentially to ensure an overlap in np of adjacent active areas. The active length in z 
is 35.6 cm for both layers. 

Data taken with the HI detector are reconstructed online and the reconstruction modules also 
determine important parameters like average drift velocities, Lorentz-angle and Tq of the CJCl and 
CJC2. In addition beam parameters like the rip position of the primary interaction point are recorded. 

In the alignment procedure, using the program Millepede, a large number of parameters is deter- 
mined; these parameter are in general corrections to detector parameters and represent small corrections. 
It turned out that a common alignment of the drift chamber and the silicon detector is essential, after an 
internal detector alignment, which already determines reasonable starting values for alignment parame- 
ters. For each of the drift chambers CJCl and CJC2 14 global parameters representing an overall shift or 
tilt are introduced. Local variations of the drift velocity VdrUt for cells halfs and layers halfs are observed, 
which are parametrized by 180-1- 112 corrections, which change with the HV configuration. For each 
wire group (8 wires) corrections to Tq are introduced (330 corrections). For the silicon vertex detector 
CST an internal alignment has already been done, using the same techniques, with in total 384 local 
parameters. In the common alignment fit 5 parameters per half ladder representing shifts are introduced 
(320 parameters). Table |l] contains all parameters. A far as possible the parameters are defined with the 
dimension length. For example the time-zero parameters are multiplied by the mean drift time. Angles 
are parametrized as shifts relative to a normalized length parameter. For example Zr is defined as z/ Zmax 
and has a range -1 to -i-l. The shift Ax with z is parametrized as Ax/zr and has the dimension length. 

In addition to physics events (ep interactions) also cosmic tracks with and without i?-field were 
simultaneously used as input to the alignment. With combined fits to both halfs of cosmic tracks parts 
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Fig. 2: The standard deviations for thie thtree 
track parameters dca (distance of closest 
aproach to the beam axis), the azimuthal 
angle ip and of the inverse transverse mo- 
mentum 1/pt are shown as a function of 
logj^Q pt/GeV, before (old) and after (new) 
the improved alignment. The standard de- 
viations are determined from the difference 
between the two parts of cosmic tracks. 



of the detectors opposite in tp angle are strongly correlated. This allows an accurate interpolation from 
the drift chamber to the silicon tracker with an accuracy a = 30 to 40 fim, which is of the same order 
of magnitude as the silicon tracker accuracy. In total about 50 000 events are used as input. Constraints 
require e.g. a zero overall shift of the detector. The CJC2 is fixed (first five parameters) vs. the cryostat 
on the basis of external survey data as reference. 

An overview about the corrections determined in the alignment fit is given in Table |l|, which has 
29 rows of different corrections. In detail these corrections are as follows. 

Rows 1-6 These rows contain geometrical parameters; in an external survey the positions are measured 

with an accuracy of a few 100 /xm and angles with an accuracy of a few 10 /irad. This accuracy is 

increased by a large factor in the alignment fit. 
Row 6 This parameter is determined from B = cosmics (straight line track model). 
Rows 11 - 14 The parameters are determined only as a check. Their actual values are determined from 

hit triplets as a function of angle /3 and Zr {(3 is the angle between the track and the normal to the 

drift direction). 

Row 15,16 The parameters change with the HV configuration. 

Row 18,20 Corrections for the wire position in and perpendicular to the drift direction. These correc- 
tions are determined in a cross check with B = Q cosmic tracks, otherwise the data from a survey 
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of flanges are used (< lO/iin). 

Row 21 These correction are used to describe the changes of drift velocity Vdnti by variations of the 
electron beam current le due to space charge effects; they are determined once. 

Row 22 - 24 These parameters change with repair/breakdown of HV cards. They are introduced in 
order to correct for bad cathode HV. 

Row 25 - 29 Corrections for the silicon vertex detector, five per half ladder, determined in and perpen- 
dicular to ladder, on average and as function of z and of ip. 

The alignment is done for run ranges and repeated after e.g. a change of the HV configuration. Run-to- 
run variations of parameters (e.g. the average drift velocity Udrift) are taken into account by the data from 
the online reconstruction as mentioned above. 



track fit 


cr(dca) [cm] 


a{l/pt) [(GeV/c)-i] 


remark 


only 1 track with Si hits 

both tracks with at least 1 Si hit 

both tracks with at least 2 Si hits 


0.0209 
0.0032 
0.0028 


5.58 • IQ-^ 
2.26 • 10-^ 
2.17 • 10^3 


drift chamber resolution 
Si tracker resolution 
Si tracker resolution 



Table 2: Standard deviations in the parameters dca and 1/pt, determined from the difference between the two parts of cosmic 
tracks. 



The alignment procedure allows to reach the intrinsic resolution not only locally but also for 
complete tracks. For the drift chambers CJCl and CJC2 the local hit resolution can be determined 
from local hit triplets as a function of drift distance; the minimum is about 120 /xm for 1 cm drift 
distance. The minimum value from track residuals is about 125 /xm. For the CST the intrinsic point 
resolution of 15 fim from hits overlaps is also reached globally. Figure ^ shows the improvement for 
track parameters, determined from data in the comparison of parameters between the two parts of cosmic 
tracks. A significant improvement is visible for larger track momenta. In Table ^ the standard deviations 
in the parameters dca and l/pt, determined from the difference between the two parts of cosmic tracks, 
are given for three different track classes. 
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